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Abstract. We present results of high-precision timing measurements of the binary millisecond pulsar PSR 
J2145— 0750. Combining 10 yrs of radio timing data obtained with the Effelsberg 100-m radio telescope and 
the Lovell 76-m radio telescope we measure a significant timing parallax of 2.0(6) mas placing the system at 
500 pc distance to the solar system. The detected secular change of the projected semi-major axis of the orbit 
X = 1.8(6) X 10~^^ It-s s~^, where x = (op sini)/c, is caused by the proper motion of the system. With this mea- 
surement we can constrain the orbital inclination angle to i < 61°, with a median likelihood value of 46° which is 
consistent with results from polarimetric studies of the pulsar magnetosphere. This constraint together with the 
non-detection of Shapiro delay rules out certain combinations of the companion mass, m^, and the inclination, i. 
For typical neutron star masses and using optical observations of the carbon/oxygen-core white dwarf we derive 
a mass range for the companion of 0.7 Mq < m2 < 1.0 Mq. We apply evolutionary white dwarf cooling models to 
revisit the cooling age of the companion. Our analysis reveals that the companion has an effective temperature 
of Tcft = 5750 ± 600 K and a cooling age of TcooI = 3.6(2) Gyr, which is roughly a factor of three lower than the 
pulsar's characteristic age of 10.4 Gyr. The cooling age implies an initial spin period of Pq = 13.0(5) ms, which 
is very close to the current period. 
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1. Introduction 

Millisecond pulsars (MSPs) have small spin periods 
(~1.5 — 30 ms), small spin-down rates (P ^ 10~^^ s s~^), 
and are believed to be recycled by the accretion of mass 
from an evolving companion star (e.g. Alpar et al. 1982). 
About ~80% of the MSP population in the Galactic plane 
are members of binary systems. Most of those move in 
nearly circular orbits around low-mass companions, form- 
ing the class of low-mass binary pulsars (LMBPs). LMBPs 
have spin periods of P ^ 10 ms, small eccentricities e ^ 
10~^, and companions of mass 0.15 Mq ^ m2 ^ 0.4 M©, 
presumably helium (He)-core white dwarfs (WDs). The 
evolutionary history of LMBPs is well understood. During 
the giant phase of the ~1 Mq companion star, mass over- 
flows its Roche lobe and is accreted onto the neutron star 
(NS). The stable mass accretion causes the pulsar to spin 
up to millisecond periods and prevents a central helium ig- 
nition of the giant, resulting in a binary system with a re- 
cycled millisecond pulsar and a low-mass He-core WD (see 
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review by Bhattacharya & van den Heuvel 1991; Phinney 
& Kulkarni 1994). 

In the last decade, pulsar surveys have revealed a new 
class of binary pulsars, the intermediate-mass binary pul- 
sars (IMBPs), with companion masses above 0.45 Mq. 
Like in the LMBP case, the IMBPs have almost circu- 
lar orbits and WD companions, but show larger spin pe- 
riods P <; 10 ms and slightly higher orbital eccentric- 
ities. As helium ignition in the giant core starts above 
a core mass of 0.45 A/q (Kippenhahn & Weigert 1990), 
the companions have to be WDs with carbon/oxygen 
(CO) or oxygen/ncon/magncsium (ONeMg) cores. PSR 
J2145— 0750 is most probably a member of the class of 
IMBPs. Discovered by Bailes et al. (1994), the 16-ms pul- 
sar is in a nearly circular orbit, with an orbital period of 
6.8 days. The system shows some of the peculiar prop- 
erties of the IMBP class: (1) The pulsar companion has 
a mass of m2 > 0.43 Mq and probably even larger than 
0.51 Mq (50% C.L.), suggesting it to be a CO-core WD. 
From optical observations of the PSR J2145— 0750 system 
and using the dispersion measure (DM) distance of 500 pc, 
Lundgren et al. (1996) concluded that the companion is 
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a WD. This is confirmed by our measurement of a timing 
parallax (see Sect. 14.1(1 showing that the DM distance is 
the true distance of the pulsar. (2) PSR J2145-0750 has 
a spin period of 16 ms, which is considerably longer than 
the average period of ~4 ms of LMBPs. (3) The surface 
magnetic field of the pulsar is quite high: i?s = 6 x 10^ G. 
(4) Finally, the spin-down age of the pulsar of ~10 Gyr is 
very high and comparable to the age of the Galaxy, raising 
interesting questions about its progenitor and the initial 
spin period. 

In contrast to LMBPs the evolution of IMBP progen- 
itor systems is not fully understood. Two models for the 
formation of heavy WDs in close orbits have been pro- 
posed: (1) The system undergoes a Common Envelope 
(CE) evolution and spiral-in phase on the Red Giant 
Branch (RGB) (Tauris 1996) or Asymptotic Giant Branch 
(AGB) (van den Heuvel 1994), similar to the evolution of 
high-mass binaries. (2) A massive companion progenitor 
looses much of its envelope on the late main-sequence or 
early RGB through highly super-Eddington mass transfer 
(Tauris et al. 2000; Taam et al. 2000). 

For the particular case of PSR J2145— 0750, van den 
Heuvel (1994) proposed an evolutionary channel with a 
CE spiral-in phase on the AGB. In this scenario the pro- 
genitor system consisted of a 1 — 6 Mq donor star, with 
an orbital period such that this star overflows its Roche 
lobe on the AGB. During the CE phase, which is similar 
to that of high-mass binary pulsar progenitors, the mass 
accretion onto the pulsar is highly unstable. This results 
in an incomplete recycling process of the pulsar and can 
explain the high values for P and Bs- One possibility to 
test these evolutionary channels is to measure the NS and 
companion masses. Here, we report on new results from 
timing measurements of the PSR J2145— 0750 binary sys- 
tem. They allow us to estimate the inclination and com- 
panion mass of the binary. From the detection of a timing 
parallax we derive an independent distance estimate of the 
pulsar. 

For a pulsar - white dwarf binary system one usually 
expects that the cooling age of the WD matches the age 
of the pulsar. This is explained by the fact that the cool- 
ing and spin-down clocks start ticking at approximately 
the same time when the companion starts to contract to a 
WD and the pulsar turns on at the end of the mass trans- 
fer phase. The cooling age of the WD companion can be 
determined from its effective temperature and its mass, 
using a WD cooling model. We use the parallax distance 
along with optical observations to revisit the cooling age 
of the WD companion, based on evolutionary WD cooling 
models. We discuss the implications of our findings for the 
evolutionary history of the binary. 



2. Observations 

We have made regular high-precision timing observations 
of PSR J2145— 0750 over a 10 yrs time span using both 
the 100-m Effelsberg radio-telescope of the Max-Planck- 



Institut fiir Radioastronomie in Bonn, Germany, and the 
76-m Lovell radio-telescope at Jodrell Bank, UK. 

The Effelsberg observations have been carried out since 
April 1994 at a centre frequency of 1410 MHz. Data 
have been acquired approximately once per month, with 
a typical observing time of 3 x 7 min. For the obser- 
vations we used a 1300-1700 MHz tunable HEMT re- 
ceiver with a system temperature of 30 K on the cold 
sky and an antenna gain of 1.5 K Jy~^. In order to mon- 
itor changes of the dispersion measure (DM) we occa- 
sionally collected data at 860 MHz, using an uncooled 
HEMT receiver with a system temperature of 65 K on 
the cold sky and a gain of 1.5 K Jy~^. During the pe- 
riod April 1994 - October 1996, the data were obtained 
using the Effelsberg Pulsar Observation System (EPOS) 
(Seiradakis et al. 1995; Kramer et al. 1997, 1998). Here, 
the two circular polarization signals were processed in a 
4 X 60 X 666 kHz filter bank combined with an incoherent 
hardware de-disperser. The detected left- and right-hand 
circular signals in each channel were de-dispersed on-line, 
added, and folded with the topocentric pulse period to 
produce a total power pulse profile of 40 MHz bandwidth. 
Since October 1996 the observations were made with the 
Effclsberg-Berkeley Pulsar Processor (EBPP), which cor- 
rects for the dispersion smearing of the signal using a co- 
herent de-dispersion technique (Hankins & Rickett 1975). 
In the total power mode the EBPP provides 32 channels 
for both polarizations with a maximum total bandwidth 
of 112 MHz, depending on DM and observing frequency 
(Backer et al. 1997). For PSR J2145-0750 a total band- 
width of 90 MHz was available at 1410 MHz. The output 
signals of each channel were fed into de-disperser boards 
for coherent on-line de-dispersion and were synchronously 
folded at the pulse period over 7-min integrations. To ob- 
tain a high signal-to-noise ratio polarization profile of PSR 
J2145-0750 (see Sect.Ol we used the EBPP in polariza- 
tion mode where the four Stokes parameters are available 
over 32 channels with a total bandwidth of 28 MHz. 

Lovell observations of PSR J2145— 0750 were made 
from October 1992 to January 1999 at centre frequen- 
cies of 606 and 1400 MHz. Both circular polarization 
signals were detected and incoherently de-dispersed in a 
2 X 64 X 0.125-MHz filter bank at 606 MHz and in a 
4 X 32 X 1.0-MHz filter bank at 1400 MHz. The signals 
were added and synchronously folded at the pulse period 
with a typical integration time of 1 — 3 min. Details of the 
observing system can be found in Gould & Lyne (1998). 

Both Effelsberg and Lovell data were time stamped 
with clock information from a local hydrogen maser clock 
and later synchronized to UTC(NIST) using the signals 
from the Global Positioning System (GPS). In order to 
calculate the pulse timc-of-arrival (TOA) synthetic tem- 
plates of the pulse profile were constructed for each back- 
end and observing frequency, and fitted to the observed 
profiles with template matching procedures (details can be 
found in Lange et al. 2001). TOA uncertainties were es- 
timated using a method described by Downs & Reichley 
(1983). Typical TOA errors for EBPP observations are 
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2.8 /us at 860 MHz and 1.9 /xs at 1410 MHz, whereas for 
Lovell observations we found 5.2 /us at 606 MHz and 4.5 /us 
at 1400 MHz. 



3. Timing analysis 

The combined TOAs, weighted by their individual uncer- 
tainties, were fitted to a spin-down model of the pulsar 
in a binary system with the software package tempo^, 
using the DE200 planetary ephemerides (Standish 1990). 
PSR J2145— 0750 lies close to the ecliptic plane, with an 
ecliptic latitude of only j3 — 5°.3, resulting in less accurate 
position and proper motion measurements if determined 
in the equatorial reference frame. In order to minimize 
covariances between the astrometric parameters we chose 
ecliptic coordinates for our timing model, performing stan- 
dard Monte-Carlo simulations to derive reliable values for 
the proper motion pix and /i^. Since the pulsar is moving 
in an almost circular orbit we applied the timing model 
for binary systems with small eccentricities (ELLl) us- 
ing the Laplace-Lagrange parameters ei and £2, as well 
as the time of ascending node Tasc (Lange et al. 2001). 
The three Keplerian parameters eccentricity e, epoch To, 
and longitude of periastron uj, are then calculated from 
the former three parameters. 

In the fitting procedure the TOA segments obtained 
with EPOS, EBPP, and Lovell data were fitted for a mu- 
tual offset accounting for different templates and TOA ref- 
erence points in the profiles. Using the full TOA data set at 
all frequencies we then determined the DM over the 10-yrs 
period. By excluding all TOAs with lines-of sight to the 
pulsar lying closer than 30° to the Sun we analysed the in- 
fiuence of free electrons in the solar system on the DM, but 
found no variations of the DM or other parameters larger 
than their Icr uncertainties. However, the multi-frequency 
timing fit resulted in a detection of a secular DM variation 
of<i(DM)/di = -2.2(4) X 10"" pccm-3 yr"!. Variations of 
pulsar DMs have been explained by a simple wedge model 
of electron density fiuctuations in the Galaxy, resulting in 
a square root dependence of d(jyM)/dt (x VDM (Backer 
et al. 1993). Recently, Hobbs et al. (2004) showed that 
this relation holds for a large group of pulsars, and de- 
rived a best-fit of d(DM)/di « 0.0002 VdM, with a scat- 
ter of one order of magnitude. For PSR J2145— 0750 our 
measurement of \d{YyM)/dt\ is consistent with the empiri- 
cal relation. By holding the values for DM and d{IiM)/dt 
fixed we then obtained the best-fit timing model for the 
astrometric, rotational and binary parameters using only 
the more accurate 1400 MHz TOAs from Effelsberg and 
Jodrell Bank. We scaled the TOA errors by an appropri- 
ate factor to achieve a uniform x^/^-dof = 1 for each TOA 
segment. 

From the best-fit model we obtained post-fit tinfing 
residuals that are displayed in Fig. ^ As seen from the 
figure, we obtained a good coverage of the 10-yrs timing 
baseline. Note that the EPOS and Lovell TOAs obtained 




Fig. 1. Post-fit timing residuals of the combined 
Effelsberg and Lovell 1400 MHz timing data as a func- 
tion of observing year. 
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Fig. 2. Distribution of the 1400 MHz post-fit timing resid- 
uals normalized by their uncertainties. 

from 1992 - 1999 have much higher uncertainties than 
the EBPP data which is due to the incoherent dispersion 
removal technique. In order to test the statistical prop- 
erties of our data we generated a histogram of the devi- 
ations of the timing residuals from the model as shown 
in Fig. 12 The post-fit residuals are consistent with the 
expected Gaussian distribution. The timing parameters 
of the best-fit model for PSR J2145-0750 are fisted in 
Table Q] Included in the Table are upper limits as well 
as derived parameters. The limits on the pulse frequency 
second derivative, i>, and the orbital period derivative, Pb, 
were found by allowing the extra parameters to vary one 
at a time in the global fit. Fitting for the derivatives of 
both the Laplace-Lagrange parameters, ei and £2, led to 
the upper limits on the eccentricity derivative, e, and the 
periastron rate of change, Co. 

4. Astrometry and intrinsic pulsar parameters 

4.1. Parallax, proper motion and transverse velocity 

Pulsar distances are estimated from their DM using mod- 
els for the electron density distribution in the Galaxy (e.g. 
Taylor & Cordes 1993, updated model: Cordes & Lazio 
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Table 1. Timing model for PSR J2145-0750 

Measured parameters'^ 



Ecliptic longitude, A 

Ecliptic latitude, /3 

Proper motion*', fi\ (masyr~^) 

Proper motion*", //^ (masyr^^) 

Right ascension'', a (J2000) 

Declination'^, S (J2000) 

Parallax, n (mas) 

Pulse frequency, v (s^^) 

Pulse frequency derivative, i/ (10^^® s^^) 

Pulse period, P (ms) 

Period derivative, P (10"^" ss"^) 

Epoch (MJD) 

Orbital period, Pb (d) 

Projected semi-major axis, x (It-s) 

First derivative oi x, x (10"^'* It-s s^^) 

ei 

£2 

Epoch of ascending node, Tasc (MJD) 

Eccentricity'*, e 

Epoch of periastron"^. To (MJD) 

Longitude of periastron'^, u; (deg) 

Dispersion measure, DM (pccm""') 

Dispersion measure derivative, d(DM)/dt (10~* 



pccm ^ yr ^) 



326? 02465688(10) 

5? 3130779(8) 

-12.37(6) 

-6.8(7) 

2l''45™50'!46726 

-07° 50' 18 '.'375 

2.0(6) 

62.2958888453343(4) 

-1.15481(8) 

16.05242365965367(14) 

2.9757(2) 

50800.0 

6.8389025099(4) 

10.1641056(6) 

1.8(6) 

-0.00000677(8) 

-0.00001805(10) 

50802.29811051(8) 

0.00001928(10) 

50806.108(4) 

200.5(2) 

9.0031(2) 

-2.2(4) 



Pulse frequency second derivative, \i>\ (s 
Orbital period derivative, |Pb| (ss^^) 
Eccentricity derivative, |e| (s^"*^) 
Periastron rate of change, \Lb\ (degyr^^) 



Measured upper limits'^ 
"^ 



< 8 X 10"^' 

< 2.4 X 10" 

< 2.8 X 10" 
<0.47 



Derived parameters 



Galactic longitude, I 

Galactic latitude, b 

Parallax distance (pc) 

Composite proper motion, /i (masyr" 

Transverse velocity, Vt (kms^^) 

Mass function, /m (Mq) 

Pulsar characteristic age (Gyr) 

Pulsar surface magnetic field (10* G) 

Timing RMS (^s) 



47? 78 
-42? 08 

14.1(4) 

33(9) 

0.024105570(4) 

10.4(5) 

6.35(14) 

2.7 



"^ Uncertainties quoted are in the last digit (s) and represent 2a estimates (twice the formal tempo errors) 
*" calculated from Monte-Carlo simulations 
'^ calculated from ecliptic coordinates 
** calculated from the ELLl model 
" Upper limits represent 95% C.L. 



2002, hereafter NE2001). Independent distance measure- 
ments define an absolute distance scale and hence are of 
fundamental importance to calibrate models like NE2001. 
For pulsars, distances are derived from trigonometric par- 
allaxes using VLBI observations, or from pulsar timing. In 
pulsar timing, the parallax is obtained by measuring the 
time delay of the TOA caused by the curvature of the radio 
wave fronts. However, the amplitude of the corresponding 
residual term, Af^ (^s) = o? cos j3 / {2cd) = 1.21 cos/?/ d, 
is very small, allowing us to detect a timing parallax only 
for a very few millisecond pulsars near the ecliptic plane. 



Here, a is the Earth-Sun distance in kpc, j3 the ecliptic 
latitude, and d the distance to the pulsar in kpc. Timing 
parallaxes have been detected for pulsars such as PSR 
B1855+09 (Kaspi et al. 1991), PSR J1713+0747 (Camilo 
et al. 1994a), PSR J0437-4715 (Sandhu et al. 1997) and 
PSR J1744-1134 (Toscano et al. 1999). Nevertheless, 
due to an ecliptic latitude of only (3 — 5°.3 we were 
able to detect for the first time a timing parallax for 
PSR J2145-0750 of tt = 2.0(6) mas. This corresponds 



to a distance estimate of d^ 



500] 



well with the DM distances, i.e. d 



120 
DM 



pc and agrees very 
— 500 pc from the 
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Taylor & Cordes (1993) model and cIbm — 570 pc from 
the NE2001 model (with typical uncertainties of ^10%). 
Using the parallax distance we derive a transverse veloc- 
ity of Vt = 33(9) kms""'^, which is in excellent agreement 
with scintillation velocities of 31 ± 25 km s^^ obtained by 
Nicastro & Johnston (1995) and Johnston et al. (1998). 

With a galactic position angle of 289° the pulsar is 
moving towards the Galactic plane. This is not in contra- 
diction to the idea that pulsars are born in the Galactic 
plane and then move to higher latitudes, as the MSP pop- 
ulation is very old (^^1 Gyr) and has already reached 
a dynamic equilibrium in the Galactic gravitational po- 
tential (e.g. Toscano et al. 1999). In this scenario PSR 
J2145— 0750, with a characteristic age of 10 Gyr (see Sect. 
I4.3|l . would have already performed many oscillations in 
the Galactic potential described by Kuijken & Gilmore 
(1989). 

4.2. Doppler corrections 

Accelerations in the Galactic gravitational potential in- 
troduce contributions to the observed period derivative 
P which cannot be neglected for MSPs that have spin- 
down rates of about six orders of magnitude lower than 
those of normal pulsars. Doppler effects arise due to (a) the 
Galactic differential rotation, (b) the vertical acceleration 
in the Galaxy (Damour & Taylor 1991), and (c) the trans- 
verse velocity of the pulsar (Shklovskii 1970). Following 
these authors, the Doppler correction to the period deriva- 
tive can be written as 




a, sin 6 



cRo 



■ cos b cos I 



f3 



sin^ I + /32 



c 



where Rq and vq denote the galactocentric radius and 
galactic circular velocity of the Sun, respectively, I and 
b the galactic longitude and latitude of the pulsar, and 
f3 — d/Ro — cos I. We choose Rq = 8.0 kpc and vn = 
220 km s^^ (Reid 1993). a^ is the vertical component of 
the galactic acceleration and can be calculated from the 
Kuijken & Gilmore (1989) model of the Galactic potential. 
For PSR J2145— 0750 we derive a kinematic correction of 
(P/P)° = 3.3(7) X 10"^^ s-\ where the major contribu- 
tion (73%) comes from the Shlovskii term {fi^ d)/c. 

4.3. Characteristic age and magnetic field 

Subtracting the kinematic contributions from the ob- 
served spin-down rate of PSR J2145— 0750 we obtain an 
intrinsic spin-down rate of P'"*'' — p — p^ — 2.45(11) x 



10 



-20 



s s which is 18% smaller than the observed value. 



It is usually assumed that the evolution of the spin 
frequency v = 1/P can be described by a power-law 
u ex —I/", where n is the so-called braking index. Then 
the pulsar age can be calculated from 



P 



(n-l)P 



'-'$ 



(2) 



where Pq is the initial period of the pulsar. Under the 
assumptions that Pq <C P and that the spin-down is 
due to magnetic dipole braking (n — 3), the spin-down 
age is given by Tq = P/(2P), the characteristic age of 
the pulsar. Subtracting the Doppler contributions, for 
PSR J2145-0750 we thus obtain t^ = P/(2P'""') = 
10.4(5) Gyr. As the uncertainty of Tc caused by the a 
priori unknown n and Pq can be considerable, an inde- 
pendent age estimate, as given by the cooling age of the 
pulsar companion (see Sect.lHll, is highly desirable. 

Assuming a pure dipolar magnetic field and using the 
intrinsic spin-down rate we can estimate the surface mag- 
netic field of the pulsar to be Ps = 3.2 x 10^ VpP'"*'' G 
= 6.35(14) X 10® G. This magnetic field is more than three 
orders of magnitude smaller than for normal pulsars, a 
clear indication than PSR J2145— 0750 has gone through 
a recycling phase. 



4.4. Timing instabilities 

Any remaining non-gaussian noise in timing residuals is 
generally explained by rotational instabilities intrinsic to 
the pulsar, or by an incomplete timing model, such as 
unmodelled DM variations or the existence of a planet in 
the system. For PSR J2145-0750 a fit for i> over the full 
timing baseline resulted in t> = —5(3) x 10"^® s~^. As this 
measurement is hardly significant we conservatively quote 
i> as an upper limit in Tabled At the same time, we do not 
detect any DM variations that are not accounted for in the 
jbjil^ng model. A simultaneous fit of our multi-frequency 
data for third- and fourth-order derivatives of DM as well 
as for i) resulted in a similar i) value, while d'^ (DM) /dt'^ 
and d'^(DM)/dt^ were not significant. 

Arzoumanian et al. (1994) analysed the timing noise 
of a number of slow and millisecond pulsars by quan- 
tifying the noisiness with the noise parameter Ag — 
log (lOl t^ / {6i/)) , where t = 10® s. They quote a fit-by- 
eye to the data of the form Ag = 6.6 -I- 0.6 log P, with a 
large scatter of the data points. Lommen (2002) revised 
this relation using new MSP data and showed that at least 
seven MSPs have significantly lower Ag's than presented 
by Arzoumanian et al. (1994). 

In order to estimate the timing noise of PSR 
J2145— 0750 we divided the total data set into three sub- 
intervals, each spanning 10® s ~ 3.2 yr of TOAs. For each 
sub-interval we transformed epoch, pulse frequency v, and 
epoch of ascending node Tasc to an epoch near the centre of 
each sub-interval. We then fitted the timing model of the 
pulsar, holding DM and d{DM)/dt fixed and allowing i) 
to vary. The resulting Ag's for each sub-interval are —4.0, 
—4.1, and —4.3, i.e. slightly larger than the predicted value 
of Ag = -5.1. We thus conclude that PSR J2145-0750 
is a possible candidate for intrinsic timing noise. Future 
timing observations will resolve this assertion. 
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4.5. Secular changes in orbital parameters 

For the first time we detected a secular variation of the 
projected semi-major axis of PSR J2145— 0750 of i = 
1.8(6) X 10"^'* It-s s~^. In order to test the robustness 
of this measurement we divided the total set into five sub- 
intervals, each spanning two years of TOAs. For every sub- 
interval the TOAs were fitted to the timing model of the 
pulsar, with the variation of x held fixed at zero. Hereby 
we transformed epoch, pulse frequency v, and epoch of 
ascending node Tasc to an epoch near the centre of each 
sub-interval. The fits resulted in an increase of x with 
time, with a slope of 1.7(9) x lO"""^"* It-s s~^, confirming 
the value from the global fit. 

The timing data showed no evidence for secular 
changes in other orbital parameters. In principle, the bi- 
nary system could undergo both a general relativistic 
(GR) advance of periastron and an orbital decay due to 
the emission of gravitational radiation (Damour & Taylor 
1991), i.e. 



^G« 



iGR 



-'£ 



1927r 



-5/3 



(TqM) 



2/3 



(3) 



P^" = —^ ^ 



(I + 73 2 , 37 4^ 



(1 - e2)7/2 



T-,5/3 i\,r-l/3 

Ip. 77117712 Ai ' 



(4) 



where the masses are in Solar units and Tq = {G Mq/c'). 
If we assume a pulsar mass of r?7i = 1.4 M© and a 
companion mass of m2 = O.SIM© (using the mass func- 
tion and a median inclination < i >= 60°), we expect 



LO 



GR 



0°.0124 yr-i and P^ 



GR, 



-5.1 X 10-16 s s" 



However, as PSR J2145— 0750 has a very small eccentric- 
ity, all effects of secular variations in oj are fully absorbed 
by the redefinition of the binary period and are therefore 
not observable in this system. The observed upper limit 
for Pb is four orders of magnitude higher which makes any 
detection of the GR orbital decay highly unlikely. 

Moreover, the Shklovskii effect discussed in Sect. 14.21 
also affects the observed value of Pb. Using the results 
from Sect. 14.21 we expect a contribution of P^'^'^ = 2.4(5) x 
10-19 s-i X P, = 1.4(3) X 10-" s s-i. This value is not 
only of different sign as the expected GR contribution but 
also about three orders of magnitude larger. Nevertheless, 
the Shklovskii term is still consistent with the observations 
as it is yet a factor of 20 smaller than our derived upper 
limited. 



5. Orbital inclination and mass 

5.1. Secular change of x 

The observed secular change of the projected semi-major 
axis X = a-p sini/c could, in principle, result from a varia- 
tion of the major axis of the pulsar orbit, Op, of the orbital 



inclination, z, or a combination of both. The contributions 
to X are (Damour & Taylor 1992; Kopeikin 1996): 



obs 



GW 



PM 



(5) 



The emission of gravitational waves (GW) leads to a 
shrinking of the pulsar orbit given by (dp/op)^^ = 2/3 x 
(Pb/Pb)^^ (e.g. Doroshenko et al. 2001). With Eq.Hone 
expects Idp/apl*^^ ~ 5.8 x 10^^^ s^^ which is seven orders 
of magnitude smaller than the observed |i:/a;|°'^^ = 1.8 x 
IQ-i^ s-i. In general, |dp/ap| will be of the same order of 
magnitude as |Pb/-Pb| for typical astrophysical processes 
in a binary. With the observed upper limit of Pb < 2.4 x 
10-12 s-i we derive |Pb/Pb| < 4.1 x IQ-i^s-i, stiU three 
orders of magnitude smaller than |i;/a;|°''''. Therefore, the 
nonzero x is most likely not caused by orbital evolution. 

The Doppler correction affects the light-travel time 
across the orbit in the same way as the pulse period, so 
that the term [x/xY-' is identical to the one for the pe- 
riod derivative (Eq. ^ see Damour & Taylor 1992). Its 
value is therefore four orders of magnitude smaller than 
the observed value {x/x)°^^. 

We conclude that the nonzero x must arise from a 
change in the observed inclination i of the orbit. The rela- 
tive motion of the binary and the observer causes a change 
of the orientation of the line-of-sight to the binary which 
results in an apparent change of i. Kopeikin (1996) derived 



PM 



= ji cot I sin i 



(6) 



where 6 denotes the difference of the position angle of 
proper motion and the position angle of ascending node. 
As is a priori unknown, we can only derive a firm upper 
limit for i using | sin^| < 1 and the lower limit of allowed 
values for x from Tabled i-e- * < tan-i(/xa;/imin) = 61°. 
We performed Monte-Carlo simulations to derive the 
distribution of inclination angles within this constraint, 
following a procedure described by Nice et al. (2001). 
Assuming that the orientation of the binary in space is ar- 
bitrary, is a uniformly distributed random variable and i 
is a random variable distributed with uniform probability 
in cosi. From these distributions we select values of 9 and 
i and retain only those combinations which satisfy Eq. El 
within the measurement errors. The resulting distribution 
for i is presented in Fig.|3 showing a broad peak towards 
higher inclinations. The median value and its asymmetric 



errors are i = 46' 



.+ 9 



^^o (68% C.L.). 



5.2. Non-detection of Shapiro delay 

The propagation of radio signals from pulsars in binary 
systems is affected by a general-relativistic time delay of 
the signals in the gravitational field of the companion 
star. For a pulsar in a circular orbit this "Shapiro delay" 
(Shapiro 1964) is given by 



At = -2 7772 Pq ln[l - sin 7 sin($ - $o) 



(7) 
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Fig. 3. Distribution of orbital inclinations i resulting from 
Monte-Carlo simulations (see Sect. l5.1|) . The median of the 
distribution is 46° (solid line), and its lower and upper la 
values are at 35° and 55°, respectively (dashed lines). 



where $ is the orbital phase in radians, $0 is the phase 
of the ascending node, and Tq — (GMq/c^). In prac- 
tice, Shapiro delay is conveniently expressed in terms of 
two observables, the "range" r — m2 Tq and the "shape" 
s = sini (Ryba & Taylor 1991), the post-Keplerian orbital 
parameters which allow a determination of the compan- 
ion mass, 7712, and the orbital inclination, i. However, for 
small inclination angles the variation of At over the orbit 
is nearly sinusoidal and cannot be separated from a small 
variation in x. For edge-on orbits {i « 90°), Ai peaks 
at $ — $0 = 7r/2, where the pulsar is behind the com- 
panion. Here, the covariance with the Keplerian param- 
eters breaks, and a measurement of Shapiro delay (and 
hence r and s) is possible. Because high inclination or- 
bits are relatively rare, Shapiro delay has been detected 
in only five pulsar-white dwarf binaries, PSR B1855-I-09 
(Ryba & Taylor 1991), PSR J1713-f0747 (Camilo et 
al. 1994a), PSR J0437-4715 (van Straten et al. 2001), 
PSR J1909-3744 (Jacoby et al. 2003), PSR J1640+2224 
(Lohmer et al. 2004), and perhaps in PSR J0751+1807 
(Nice et al. 2003). 

We did not detect Shapiro delay in the timing data of 
PSR J2145— 0750. However, we can use the non-detection 
to place limits on allowed values for i and 77^2. Following a 
procedure described by Nice et al. (2001), we tested a 2D 
grid of timing models, that include the Shapiro parameters 
in the range of < sinz < 1 and < 7712 < 1.4 Mq, allow- 
ing all other parameters to vary. The resulting x^ values 
are mostly very similar to a model without Shapiro delay. 
However, for models with high inclination angles and/or 
high companion masses we obtain significantly higher val- 
ues for x^- The failure of these models means that one 
would have detected Shapiro delay had this been the true 
values of 7 and 7772- In Fig.^ the parameter space excluded 
by models with Ax^ > 9 (i.e., 3cr) is shown. 



5.3. Polarimetry 

During the recycling process of MSPs the neutron star is 
spun-up by mass transfer from the companion star which 



is expected to cause an alignment of pulsar spin and or- 
bital angular momentum of the binary system (Phinney 
& Kulkarni 1994). Assuming that this holds for PSR 
J2145— 0750 it is then possible to constrain the orbital 
inclination from the known geometry of the pulsar mag- 
netosphere. We therefore observed PSR J2145— 0750 at 
1410 MHz using the EBPP in polarization mode and de- 
rived a position angle that is similar to the findings of 
Sallmen (1998) and Stairs et al. (1999). The profile of 
PSR J2145— 0750 shows at least three components, with a 
bridge of emission between the two strongest components. 
Unfortunately, the linear intensity at 1410 MHz is very 
weak (but see also Xilouris et al. 1998; Sallmen 1998), 
and the accompanying position-angle swing is extremely 
complex, with multiple orthogonal mode changes (Stairs 
et al. 1999). Applying a rotating vector model (RVM) 
(Radhakrishnan & Cooke 1969) to the observed position 
angle usually leads to high uncertainties in the RVM fit. 
We fitted a RVM curve to the position angle of the whole 
profile (including the precursor) and found a best fit for 
a magnetic inclination of a ~ 52° and an impact angle 
of cr ~ —3°. These values are considerably different from 
those found by Sallmen (1998) which is probably due to 
the inclusion of the precursor in our fit leading to a longer 
baseline in pulse longitude. Even though the position an- 
gle is poorly constrained leading to considerable errors in 
the RVM fits, applying our best-fit values the angle be- 
tween the pulsar rotation axis and the line-of-sight results 
in a + a ^ 49°. If one assumes that the pulsar rotation 
axis is aligned with the orbital momentum, one derives 
i = a + a{ori = 180° - {a + a)), i.e. i ~ 49°. 

One generally needs to be careful in interpreting the 
position angles observed in recycled pulsars as being the 
reflection of the underlying geometry. The position angle 
swings are typically much flatter and often do not resem- 
ble the expected S-like swing (e.g. Xilouris et al. 1998; 
Stairs et al. 1999). Given this caveat in mind, it is inter- 
esting to note that the fitted inclination angle agrees very 
well with the median value of i = 46° obtained from our 
X measurement. We conclude that the analysis of the pul- 
sar magnetosphere by polarimetry results in a consistent 
picture of the orbital inclination, strongly supporting our 
earlier findings. 



5.4. Companion mass 

As demonstrated, the x measurement and the limits on 
Shapiro delay for PSR J2145— 0750 can be used to put 
constraints on the orbital inclination, i, and the compan- 
ion mass, m,2, as shown in Fig.^ From the observed x we 
derive an upper limit on the inclination angle of 7 < 61°, 
so that all points with cos 7 < 0.48 are excluded as indi- 
cated by the hatched rectangle in the figure. The curved 
hatched area is the parameter space excluded due to the 
non-detection of Shapiro delay in the binary. From the 
mass function of PSR J2145— 0750 we derive the displayed 
contours for constant neutron star masses. 
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Fig. 4. Constraints on orbital inclination, i, and WD mass, m2, from i-measurement and limits on Shapiro delay. 
Points inside the hatched rectangle are excluded by measured i, points inside the curved hatched area are excluded 
by the non-detection of Shapiro delay. From the mass function we calculated the contours for constant neutron star 
masses, indicated by the solid curves. The vertical lines denote the median value {i = 46°, solid) and its la errors (at 
35° and 55°, dashed) for the distribution of inclination angles derived from the observed x. From these uncertainties 
we derive the companion mass limits 0.5 Mq < m2 < 1.0 Mq displayed by the horizontal dashed lines. Photometric 
studies of the WD companion (see Sect.|SJ) result in a lower mass limit of 0.7Mq (horizontal dotted line). 



We cannot put any constraints on the companion mass 
from evolutionary considerations of the binary as the PSR 
J2145— 0750 system most probably went through a phase 
of highly unstable mass accretion (van den Heuvel 1994). 
Thus, an unique relation between the orbital period, Ph, 
and the secondary mass, 1112, as found for LMBPs, is not 
applicable for this system. However, from the distribution 
of allowed inclination angles (see Sect. l5.l)l . and the con- 
tours for constant pulsar masses we can estimate a range of 
allowed companion masses. Using recent results of neutron 
star mass measurements (Thorsett & Chakrabarty 1999; 
Stairs 2004 and references therein) we apply a range of 
1.2 Mq < nil < 1.8 A/q for the pulsar mass (for a discus- 
sion see Sect.O. As shown in Fig.01 for inclination angles 
within their la uncertainties we derive a range of allowed 
companion masses of 0.5 Mq < TO2 < 1.0 Mq. 

6. The nature and cooling age of the white dwarf 
companion 

From the timing analysis we derive a range of companion 
masses of 0.5 Mq < 1112 < 1.0 Mq, and a global mini- 
mum of 7712, min — 0.48 Mq for an extremely low pulsar 
mass of nil — 1.2 Mq (see Sect. [Tj). The upper mass 
limit for a helium (He)-core WD with solar metallicity 



is MHc-WD,max ~ 0.45 Mq. Hcuce we conclude that the 
companion is a WD with carbon/oxygen (CO) core in 
the intermediate WD mass range. This result is consistent 
with a CE evolution of the binary system as proposed by 
van den Heuvel (1994). Using the derived mass range of 
the WD companion and its photometrically determined 
effective temperature it is possible to determine an age 
of the binary system. This age is based on WD cooling 
models coming from stellar evolution theory and can be 
compared with the characteristic pulsar age. 

Optical observations of the WD companion with the 
HST from Lundgren et al. (1996) yielded my = 23.7 ±0.1 
and 77ii = 23.0±0.1. From these measurements the authors 
derived an effective temperature of T^s — 5800 ± 300 K 
using the color-temperature calibration from Bergeron et 
al. (1995) and the bolometric correction from Monet et al. 
(1992). With the same observational data from Lundgren 
et al. (1996), Hansen & Phinney (1998) estimated T^s = 
6235 ± 970 K based on the V — I broad-band colors of 
Bergeron et al. (1995). As conservative limits, we thus 
use a range of effective temperatures for the companion 
between 5250 and 7200 K. 

Fig. O shows an age - Toff diagram with evolution- 
ary tracks of CO-core WD models from Blocker (1995) 
and Blocker et al. (1997). The horizontal lines indicate 
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Fig. 5. Age - temperature diagram with cooling curves 
from evolutionary models of Blocker (1995) and Blocker 
et al. (1997) for CO-core WDs (0.52 Mq < Mwd < 
0.94 Mq). The labels give the WD masses in M©, and 
the straight horizontal lines indicate the possible range 
of effective temperatures for the WD companion of PSR 
J2145— 0750 according to the studies of Lundgren et al. 

(1996) and Hansen & Phinney (1998), i.e. 5250 K < Tcs < 
7200 K. The thick dashed vertical lines mark the intersec- 
tions of these lines with the cooling curves from the least 
and the most massive CO-core model (AfwD = 0.524 and 
0.940 Mq, respectively) to indicate the corresponding age 
range. 

the observational range for Toff- The curve at the very 
left of the diagram belongs to the CO-core WD model 
with the lowest mass, followed by curves of increasing core 
mass while going to right. The WD models with masses 
of 0.524 and 0.940 Mq roughly bracket the mass range 
for the companion derived from the timing analysis. As 
the vertical dashed lines in Fig. illustrate, a WD with 
mwD = 0.524 Mq would have an age of 1.3... 2. 6 Gyr, 
whereas a massive WD with mwo — 0.940 Mq would have 
an age of 2.8 .. . 4.9 Gyr. 

Using the photometric measurements of Lundgren et 
al. (1996) in the B, V, and / bands together with our par- 
allax distance for PSR J2145— 0750, the mass, age, and 
effective temperature of the WD companion can be fur- 
ther constrained. With d = 500ti2o P'^ ^^'^ ^^^ appar- 
ent magnitudes tob = 23.9, my = 23.7, and nii — 23.0 
from Lundgren et al. (1996) we find absolute magnitudes 
Mb = 15.4l^j, My = 15.2l°j, and Mj = 14.5l^j. 
Applying the bolometric corrections BCi — 0.5 from 
Monet et al. (1992) and BCy = -0.2 from Bergeron et 
al. (1995), from both My and Afj we obtain a bolometric 
brightness of Afboi = 15.0^Qg, which finally gives a WD 
luminosity of log(L/LQ) — — 4.1^q'2. This value can be 
compared with the luminosity predictions from the WD 
cooling models from Blocker (1995) and Blocker et al. 

(1997) for Tcff == 5250 . . . 7200 K. The model parameters 
are summarized in Table |21 

As can be seen from Tab.EJ using log(i/LQ) — —4.1 as 
an additional constraint a CO-core WD with m,2 ^ 0.7 Mq 



Table 2. Cooling properties of the CO-core WD 
companion of PSR J2145— 0750. All values are taken 
from WD evolutionary models of Blocker (1995) and 
Blocker et al. (1997) 



Model parameters of the WD companion 
Mwd /Me \og{L/LQ) Teg/K r,ooi7G^ 



0.524 
0.524 



-3.348 
-3.913 



7200 
5250 



1.35 
2.60 



0.605 
0.605 



-3.448 
-3.995 



7200 
5250 



1.50 
3.00 



0.696 
0.696 



-3.532 
-4.087 



7200 
5250 



1.90 
3.70 



0.836 
0.836 



-3.661 

-4.285 



7200 
5250 



2.40 
3.40 



0.940 
0.940 



-3.825 
-4.850 



7200 
5250 



2.80 
4.85 



can be ruled out as companion for PSR J2145— 0750, since 
WDs with lower masses would be too luminous. Here, a 
WD withm2 « 0.7 Mq would just havelog(L/LQ) = -4.1 
in accordance with the observations, and the cooling age 
would be Tcooi = 3.7 Gyr. On the other hand, a WD with 
7712 = 0.94 Mq would be too faint at a temperature of 
Tcff — 5250 K. Using log(L/LQ) — —4.1 as constraint, we 
obtain TcooI = 3.4 Gyr from the 0.94 Mq cooling model, 
which in turn leads to a new upper limit for the effec- 
tive temperature of Tes = 6400 K (see Fig.|5ll. Therefore, 
from the overall comparison of observational constraints 
and WD models, we conclude that the PSR J2145-0750 
companion has a mass of 0.7 Mq < m2 < 1.0 Mq, an ef- 
fective temperature of T^g = 5750 ± 600 K, and a cooling 
age of Tcooi = 3.6(2) Gyr. 

The WD mass range derived here is in accordance with 
the results from our timing analysis (see Sect. 15.4(1 . but 
provides a tighter lower limit of 0.7 Mq. On the other 
hand, the results on the WD properties are also in agree- 
ment with previous studies from Lundgren et al. (1996) 
and Hansen & Phinney (1998). Taking the mass range 
0.7 AfQ < m2 < 0.9 Mq, we can also compare our re- 
sults with those obtained from the recent WD models 
of Richer et al. (2000). Using their models, we obtain 
Teff = 5700 ± 350 K and TcooI = 4.4(1) Gyr in line with 
our results. 



7. Discussion &. Conclusions 

Our timing observations of the PSR J2145— 0750 binary 
system constrain the orbital inclination angle to be i < 
61°, with a median value of i = 46°. Polarimetric stud- 
ies of the pulsar magnetosphere lead to a consistent result 
but cannot provide stronger constraints on the orbital in- 
clination. From a statistical analysis of allowed inclina- 
tion angles we are able to derive proper mass limits for 
the pulsar companion if we apply a realistic estimate for 
the pulsar mass. Neutron stars are believed to have mass 
close to the Chandrasekhar value of ~1.35Mo. Thorsett 
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& Chakrabarty (1999) reviewed all neutron star mass 
measurements and found a remarkably narrow distribu- 
tion of 1.35(4) Mq. Recent observations revealed neutron 
stars that have lower masses, e.g. PSR J0737— 3039B with 
1.250(5) A/q in the newly detected double-pulsar system 
(Lyne et al. 2004), and PSR J1141-6545 with 1.30(2) Mq 
in a NS-WD system (Bailes et al. 2003). New timing obser- 
vations showed that recycled pulsars in NS-WD binaries 
tend to have masses that are significantly higher than the 
canonical value of 1.35 Mq (Nice 2003), which is generally 
explained by a long phase of extended mass transfer from 
the companion to the NS. This also holds true for PSR 
J0621-H1002 with a pulsar mass of ~1.7M0 (Splaver et 
al. 2002), an IMBP system with properties similar to PSR 
J2145— 0750. The currently known NS masses in recycled 
pulsar- WD systems lie in the range 1.5 — 1.7 Mq, with 
a median value of --I.6M0 (ignoring PSR J0751-I-1807 
due to the large uncertainties in its mass measurements; 
Stairs 2004 and references therein). Taking into account 
the considerable uncertainties, we conservatively choose a 
range of 1.2 Mq < mi < 1.8 M© for the pulsar mass of 
PSR J2145-0750. As discussed in Sect. El this results in 
allowed companion masses of 0.5 AIq < 7712 < 1.0 AIq. We 
hence conclude that the pulsar companion is a CO-core 
WD in the intermediate WD mass range. For the median 
pulsar mass of ^1.6 Mq and the observed median inclina- 
tion angle of 46° we obtain a median companion mass of 
~0.7 Mq. 

Using the measured parallax distance and the con- 
straints coming from photometrical observations of the 
WD companion of PSR J2145-0750, the WD proper- 
ties can be further constrained. Our analysis reveals that 
the companion has a mass in the range 0.7 Mq < m,2 < 
1.0 Mq, a cooling age of TcooI = 3.6(2) Gyr, and an effec- 
tive temperature of Tcff — 5750±600 K, in good agreement 
with results from previous studies (Lundgren et al. 1996, 
Hansen & Phinney 1998, Richer et al. 2000). In particu- 
lar, the WD model for a companion mass of 0.7 Mq, the 
median value derived from the timing analysis, exactly re- 
produces the measured WD luminosity for Toff = 5250 K. 

The WD cooling age is roughly a factor of 3 lower than 
the pulsar's characteristic age of 10.4 Gyr, which is gen- 
erally considered as a rough upper limit of the system age 
(Camilo et al. 1994b). Under the assumption that the pul- 
sar's magnetic field does not decay, we can use the WD 
cooling age to estimate the initial spin period Pq of PSR 
J2145-0750, by inverting Eq. O The Doppler contribu- 
tions are subtracted as we apply the intrinsic spin-down 
rate /""*"■. For the derived system age of r = twd = 3.6(2) 
Gyr and typical braking indices n = 2. ..3 (Lyne 1996 and 
references herein) we hence obtain an initial spin period 
of Pq = 13.0(5) ms, which is very close to the current pe- 
riod. Theories regarding the spin-up of recycled pulsars 
to millisecond periods (Smarr & Blandford 1976; Ghosh 
1995 and references therein) predict that the initial pe- 
riod should be equal to the equilibrium spin period of 
the neutron star in the X-ray binary phase, accreting at 
the Eddington rate MEdd- For PSR J2145-0750 we find 



M/MEdd = 0.004(1). Hence the final accretion rate of the 
pulsar was highly sub-Eddington, which was also observed 
in other IMBP systems (Lundgren et al. 1996). This is not 
surprising as the unstable CE phase of the binary was too 
short to spin up the pulsar to the equilibrium spin period. 
Instead, as proposed by van den Heuvel (1994), the initial 
period may have resulted from sub-Eddington accretion 
from the stellar wind of the giant star before reaching the 
asymptotic giant branch phase. 

The prospects for measuring the pulsar mass of the 
PSR J2145— 0750 system are quite low. As the pulsar com- 
panion is bright enough in the optical, the system masses 
could in principle be determined by a spectroscopic study 
(see van Kerkwijk et al. 1996). However, given the com- 
position of the companion, suitable lines might not be 
present. Similarly, post-Keplerian effects beyond Shapiro 
delay will not be detectable in the foreseeable future. As 
mentioned in Sect. 14.51 relativistic effects such as orbital 
period decay or periastron advance are fully absorbed by 
the redefinition of the binary period and therefore not ob- 
servable in this low-eccentricity binary. The precision of 
the Shapiro delay measurement (or limit) cannot be im- 
proved drastically with existing radio telescopes. It will 
need telescopes like the Square-Kilometre- Array (SKA) to 
provide a large sample of precision mass measurements. 
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